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Abstract 

In high energy heavy ion collisions a hot and dense medium is formed, where the Ua(1) or 
chiral symmetry may temporarily be restored. As a consequence, the mass of the ?/(958) mesons 
may be reduced to its quark model value, and the abundance of rf mesons at low px may be 
enhanced by more than a factor of 10. The intercept parameter A* of the charged pion Bose- 
Einstein correlations provides a sensitive observable of the possibly enhanced rf abundance. We 
have analyzed A*(?t7.t) data from ^s NN = 200 GeV central Au+Au reactions measured at the 
BNL Relativistic Heavy Ion Collider (RHIC), using extensive Monte Carlo simulations based on 
six popular models for hadronic multiplicities. Based on the combined STAR and PHENIX data 
set, and on various systematic investigations of resonance multiplicities and model parameters, we 
conclude that in y/s NN = 200 GeV central Au+Au reactions the mass of the rf meson is reduced by 
Am*/ > 200 MeV, at the 99.9% confidence level in the considered model class. Such a significant rj 
mass modification may indicate the restoration of the V a(1) symmetry in a hot and dense hadronic 
matter and the return of the ninth "prodigal" Goldstone boson. A similar analysis of NA44 S+Pb 
data at top CERN Super Proton Synchrotron (SPS) energies showed no significant in-medium rf 
mass modification. 
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I. INTRODUCTION 



In terms of the quark model, one can observe a spontaneous symmetry breaking of the 
approximate SU(3)l x SU(3)r symmetry, resulting in nine pseudo-Goldstone bosons, that 
are usually associated with the light mesons formed as u, d, s quark- ant iquark bound states. 
This naive picture is, however, complicated by the fact that the rj' meson has a large mass of 
the order of 1 GeV. As early as 1970, Kobayashi and Maskawa concluded that the large mass 
of the rj meson (formerly known as X) is a serious problem that is difficult to understand in a 
chiral SU(3)l x SU(3)r model with an explicit symmetry breaking term between singlet and 
octet states [1]. They found that the existence of an effective six-quark determinantial vertex 
is necessary. As shown by 't Hooft in 1976, this vertex is contained in instanton-induced 
quark interactions [2]. An interesting aspect of this Kobayashi-Maskawa-'t Hooft or KMT 
term [3] is that it can give rise to a flavor mixing in the scalar as well as in the pseudoscalar 
channels. The coupling between the pseudoscalar singlet and octet states r)o and r/g arises 
both from the SU(3)y breaking and the anomaly terms, assuming isospin symmetry. The 
physical r\ and rj mesons are given by the mixing of the rjs and 7? modes, and the mass of 
the rjo singlet state turns out to be sensitive to the strength of the KMT vertex. An explicit 
calculation for the general case gives the mixing angle 9 v (mi) = —20.9° [3]. 

The rj and rj mesons change their masses as a function of the temperature T, due to 
both the T dependence of the quark condensate, and the possible decrease of the KMT 
coupling constant with increasing T. The mixing angle 8 V also becomes T dependent: as 
the temperature increases, mixing between 77 and rj approximates the ideal one, and the 
r/o component in the physical 77' decreases. On the other hand, with increasing T, the 770 
tends to play the role of the ninth Nambu-Goldstone boson of the SU(3)l x SU(3)r x C/a(1) 
symmetry, and loses its mass rapidly. At low temperatures, the ?7a(1) P ar t of the symmetry 
is broken by instantons, invoking distinct vacuum states. Tunneling between these vacuum 
states is only possible "at a cost", giving extra mass to the 77' meson. However, as the 
transition amplitude is dependent on the strong coupling constant a s , it follows that the 
effect of instantons rapidly decreases with increasing energy density. This is an effective 
restoration of the Ua(1) symmetry at finite T, first suggested in Ref. [4]. Thus in high 
energy heavy ion collisions, where a hot and dense medium is created, the Ua(X) symmetry 
may temporarily be restored [5-7]. 
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In Refs. [8, 9] we reported the first observation of a significant reduction of the 77' mass, 
based on an analysis of PHENIX and STAR data [10, 11] from y/s NN =200 GeV central 
Au+Au collisions at the BNL Relativistic Heavy Ion Collider (RHIC). The subject of the 
present manuscript is to detail and systematically explore these experimental signatures of 
a partial Ua(1) symmetry restoration at RHIC. 

In thermal models, the production cross sections of the light mesons are exponentially 
suppressed by the mass. Since the 7/ mesons are heavy, by default one expects the number 
of 7/ mesons to be about two orders of magnitude less than the number of pions. However, 
as a consequence of the mass reduction, this suppression would be moderated, and the 7/ 
mesons would show up in an enhanced number. Once produced, the 7/ is expected to be 
decoupled from other hadronic matter, since its annihilation cross section is very small. At 
the same time, the low-px fj mesons are trapped in the medium due to energy conservation 
reasons, forming a condensate. Due to expansion, the medium cools down and freezes out, 
emitting asymptotic particles that propagate in T = vacuum. By this time, the 7/ mesons 
regain their original mass, hence the enhancement will mostly appear at low [5-7]. It 
is to be noted that the 7/ lifetime is much longer than the lifetime of the hot and dense 
medium, therefore a direct observation of the mass shift seems to be extremely difficult. 

A promising channel of observation is the dileptonic decay 7/' — > £ + £~^, because a low-pp 
rj enhancement would give extra lepton pairs to the low invariant mass region. The paper 
of Kapusta, Kharzeev and McLerran on the return of the prodigal Goldstone boson [6] was 
in fact motivated by the dilepton enhancement seen in CERN Super Proton Synchrotron 
(SPS) Eiat, = 200A GeV/c laboratory energies in S+Pb reactions. Recent interpretations 
of CERES [12] and NA60 data [13] indicate that the approach to a chiral symmetry re- 
stored state could proceed through resonance broadening and eventually subsequent melt- 
ing, rather than by dropping masses or mass dependency or mass degeneracy between chiral 
partners [14]. Recent PHENIX findings also show a definite excess in the m e + c - < 1 GeV 
dielectron invariant mass region in y/s NN = 200 GeV Au+Au collisions [15]. Unlike at lower 
beam energies, in this case the contribution from a hot hadronic phase without mass shifts 
seems to be insufficient to account for the enhancement seen in the data [16]. 

In the present work we search for the 77' mass modification and the related restoration of 
the chiral symmetry in a different channel, using already published like-sign Bose-Einstein 
correlation (BEC) measurements as proposed by Ref. [17]. We report on a detailed system- 
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atic study of the presently available published data sets. Future high precision data points 
in the 1ow-£>t region will help to reduce the uncertainties of the measurements, and hopefully 
will provide a more precise estimation of the rj mass reduction, compared to our current 
analysis. 



II. BOSE EINSTEIN CORRELATIONS 



Correlations between pions carry important information about the space-time structure of 
the medium created in heavy ion collisions. The widths of r\ and 7/ are = 1.30 ±0.07 keV 
and Trf = 204 ± 15 keV, corresponding to large decay times: they produce pions at cr^ = 
967 fm and CT n = 152000 fm, which are huge compared to the characteristic HBT radii 
of 4-6 fm. Among the decay channels of the rj', the rj' — > r]7T + 7i~ channel has the largest 
branching ratio of about 45%. Furthermore, the r\ mesons decay into charged pions: The 
r\ — > 7T + 7r~7r° and rj — > 7r + 7r _ 7 processes together have a branching ratio of approximately 
27% [18]. 

The Bose-Einstein correlation function of pion-pion pairs can be expressed in terms of 
the relative and the mean four- momenta, Ak = A/c M = p^ — p% and K = = (p± +p%)/2, 
respectively: 

C(M,K)= ^ 2(Pl ' P2 \, (1) 
1 ; Ni(pi)JVi(p 2 ) 1 ; 

where pi^ are the three-momenta of particles 1 and 2, Ni and N 2 are the one- and two- 
particle invariant momentum distribution functions. 1 

It has been shown that the source can be handled in the core-halo picture, where the 
"core" consists of primordially created pions and those ones coming from fast-decaying 
resonances, while the other pions, coming from more slowly decaying resonances, make the 
"halo" [19, 20]. The core region is resolvable by BEC, while the halo is not. 2 In fact, the 
core— halo separation always depends on the experimental two-track resolution. For example 
in PHENIX and STAR, two tracks are separable with a momentum difference of 5Q larger 
than 4 to 5 MeV, corresponding to a spatial separability 5x that dies off at about 40-50 
fm [22, 23]. Long tails extending to this region were recently observed by the PHENIX and 
NA49 collaborations at RHIC and CERN SPS energies [22, 24]. Note that these long tails 

1 The usual m 2 — p 2 = p^p^ = E 2 — p 2 notation is used here, where p 11 = (E,p) and p is the three- 
momentum. 

2 This has been tested numerically in [21]. 



are also seen in kaon imaging, while the bulk production is well described with characteristic 
scales of 4-6 fm [25]. 

In the core-halo picture [19, 20] the correlation function can be measured with the so 
called extrapolated intercept parameter A* as 

C(Ak, K) = 1 + A*i? c , (2) 

where R c is defined by the Fourier transform of the one-pion emission function S c (x, K) of 
the core as 

_ \S c (Ak,K)\ 2 
c ~ \S c (0,K)\* ' 

S c (Ak, K) = J d 4 xS c (x, K)e iAkx , (3) 

and the intercept parameter A* is derived from the extrapolation of the correlation value 
to AA; = 0, as A* = C(0,K) — 1. However, this extrapolation does not include the corre- 
lation between halo-halo and core-halo particle pairs, supposed to be unresolvable by our 
detectors. 3 

The intercept parameter can also be expressed directly with the fraction of the core pions 
to the total number of pions: 

where 

NZa l0 (mr) = N^imr) + N^ n+ (m T ) + A^ + (m T ) + ^L + (%) , (5) 



and it is per se sensitive to the core vs. halo ratio. (In the above formulae = m 2 + p\ 
denotes the transverse mass of the 7r + s, and the A^ 7r+ (mx)'s are the mi distributions of 
the corresponding decays.) If one does a one-dimensional investigation characterized by the 



invariant momentum difference q = \J (k\ — A^) 2 , and assuming a Gaussian approximation 
of the source, one will have Eq. (2) in the form 

C(q) = 1 + \e~\ qR \ 2 . (6) 

In this case A = A*. 



For fully thermal particle emitting sources, the exact value of intercept is \ xc t = 1. Generally, \ X ct > A*. 
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As A* is defined in Eq. (2) as an extrapolation of the experimental data to Ak^ — > 0, it is 
essential that it does not depend strongly on the method of the extrapolation. The Gaussian 
assumption for the shape of the source is a model hypothesis that has been shown not to be 
precise [22]. Instead, one can choose an experimental procedure independent of theoretical 
assumptions, and get a description of arbitrary precision with the Edgeworth expansion by 
completing the Gauss shape with an infinite series in the space of Hermite polynomials [26]. 
In this case, 



C(q) = N 



I tAJ 

l + X E e-^[l + J2^H n (V2 q R E 



n=3 



(7) 



where M is an overall normalization parameter, the n n is a weight factor to the nth Hermite 
polynomial 

^(*)=W4Y>/ 2 , (8 ) 



dt / 

and Re and Xe are the Edgeworth radius and intercept parameters. In practice, an expansion 
up to the sixth order is sufficient [11]. Note that, as the even order Hermite polynomials 
have non- vanishing values at t = 0, Xe is different from A* here: 



A* — A; 



~8~ ~ 48 



(9) 



A straightforward generalization of these approximations to the three (out, side, long) 
dimension case is given in [27]. 

In the following we will show that the phenomenon of a partial {7a(1) restoration and a 
connected mass reduction of the if meson is able to explain the behavior of A* measured in 
RHIC experiments. A previous analysis, based on the same idea, was carried out on NA44 
data, and showed that a reduction of the mass should result in a dip of the A*(mx) 
at low-mx values. However, no significant signal of such a mass modification was seen in 
E lab = 200A GeV S+Pb collisions at CERN SPS (Fig. 1 of [17]). 



III. DATA SETS 



In contrast to the SPS data, a low-mx dip of the A,(™t) was measured in ^s^n — 200 GeV 
Au+Au collisions at RHIC by both the STAR and the PHENIX experiments [10, 11]. The 
A^mx) points were extracted from the measured correlation function of like-sign pions, using 
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different methods. A comparison of the measurements to the FRITIOF [28] calculation 
including variation of the 7/ mass was presented in Fig. 4. of Ref. [29]. It is important to 
note that, although the Gaussian fit typically has 1-2% error on the A value, this does not 
reflect the real error on the intercept caused by the extrapolation. For example, exponential 
fits yield larger values for A than the Gaussian fits do, and the difference between Xoauss 
and Xexponentiai is larger than several (sometimes more than 5) standard deviations [30]. A 
reasonable range of errors can be estimated with the help of an Edgeworth fit [29]. In the 
detailed presentation of Ref. [29] it was suggested to utilize the normalized \*{m^) / X^ ax 
quantity in order to remove sensitivity to the extrapolation technique, and to reduce other 
systematic errors. Note that \™ ax is defined as the value of the extrapolated intercept 
parameter in an mx region where its value is saturated; according to simulations as well as 
the presently used PHENIX and STAR measurements, this corresponds to 0.5 GeV < < 
0.7 GeV region. 4 One should also note that the A*(mx) data may depend on the goodness of 
the particle identification too: other particles misidentified as pions will reduce correlation 
and will push the measured A*(mx) data down. 

The PHENIX A*(mx) data set was derived using the Bowler-Sinyukov method [31], while 
the STAR data set was obtained from a sixth order Edgeworth fit. Here we apply the 
method of Ref. [26] using the values and errors of the Gaussian A and the K i ri Edgeworth fit 
parameters taken from [11] in order to compute the sixth order A* values using Eq. (9). We 
use complete error propagation in order to determine the corresponding uncertanities of A*. 
Although the Gaussian A and the Edgeworth A# parameters are significantly different, the 
extrapolated intercept parameter A* is similar to A within errors (see Table V and Fig. 12 
in Appendix A). However, the error bars on A* are significantly larger than the naive errors 
on A. As emphasized before, the Edgeworth A* gives a more realistic estimation of both the 
value and the error of the intercept parameter, hence of the core-halo ratio. 

Another data set on A* from STAR was obtained with a Gaussian fit using the Bowler- 
Sinyukov method, and it shows a good agreement with the STAR Edgeworth data (Fig. 14 
of [11]). We do not use this data set in our analysis, since the error bars do not include all 
the relevant systematic effects, being in the order of 5X ~ 0.001, which is way smaller than 
the systematic error coming from the choice of the extrapolation. PHENIX preliminary 
4 Here \™ ax is the A*(mx) value taken at mx = 0.7 GeV, with the exception of the STAR data, where the 

data point at the highest mx = 0.55 GeV is considered. Note that the mx dependency of the A*(mx) 

measurements in the 0.5-0.7 GeV region is very weak. 
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data [29] are also shown, for comparison purposes only, as their systematic errors are not 
yet finalized. The data sets for A* are detailed in Fig. 1 before and after normalizing with 
\™ ax . Each of these data sets indicate the dip of A* in the 1ow-pt region. 

Given that we are interested in obtaining final errors that include all relevant systematic 
effects, we decided to analyze simultaneously the PHENIX final Sinyukov-corrected A* data 
set together with the STAR Edgeworth A* data set, both normalized to their \™ ax values. 
One can note that these two data sets are not in perfect overlap with each other, although 
they are consistent within their errors, and it is possible to fit both data sets simultaneously 
with good confidence levels. The Gaussian A(mx) data set of STAR has errors that appar- 
ently do not include systematics from the Gaussian ansatz, thus a quantitative comparison 
to our model was not reasonable. However, the best fits to the other data sets qualitatively 
agree with the Gaussian A(rax) points. We have also checked that separate analysis of the 
STAR and PHENIX data sets yields results which are consistent with the presented results, 
as part of the systematic studies. However, the combined PHENIX and STAR data set 
provided a more precise estimate for the allowed regions of the model parameters. 

The difference between the selected PHENIX and STAR data sets possibly reflects the 
systematic error from different experimental conditions, for example, PHENIX data were 
measured in the 0-30% centrality class, while STAR Edgeworth results were published for 
the 0-5% centrality selection, and the particle identification in the two experiments is also 
different. These differences resulted in a systematic uncertainty of our analysis, too, however, 
as we shall demonstrate at the end, these systematic errors are still of the order of the 
statistical uncertainties and the dominant error in estimation of the in-medium modified 
mass of the rj mesons comes from the choice of the resonance model and its parameters. The 
relative systematic error from the difference in centrality selection is estimated in Appendix B 
to be not larger than 9.8%. 

IV. MODELING AND SIMULATION 

First simulations of the {7a(1) restoration at SPS S+Pb collisions at Ei a b = 200/1 GeV 
bombarding energy had predicted a dip of A*(mx) at mx values below 0.25 GeV [17]. It was 
found that the depth of this dip was governed by the value of m*, as an input parameter 
for the simulations. In those simulations, the r( mesons from the condensate had been 
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FIG. 1. (Color online) Data sets of A*(mx) (left) and X^m^) / \™ ax (right) from RHIC ^snn = 
200 GeV like-sign pion correlation measurements. Note that the errors on the STAR Gaussian 
data set (Fig. 14 of Ref. [11]) are in the order of 0.001 MeV (not all systematics included), and one 
of the PHENIX data sets has only preliminary errors [29]. Hence we are left with the remaining 
two data sets when evaluating the systematic errors on m*,: The published PHENIX data set of 
Ref. [10], and the STAR data set calculated from the data shown in Fig. 13 of Ref. [11] using the 
method of Appendix A. 



assumed to have no transverse momenta at all, resulting in a very steep hole-like structure 
of the Iow-m T part of the A*(m T ) [17]. Such an oversimplification is not adequate for the 
description of the dip in the RHIC data. To improve on it, we introduce here an effective 
thermal spectrum for the mesons from the condensate, characterized by an inverse slope 
parameter B~ x . As we demonstrate below, this parameter B~ l controls the steepness of the 
dip of A* (m^), while ml, controls its depth. Also, instead of relying on a given model of 
resonance production (as on FRITIOF [28] in Ref. [17]), here we utilize six different models, 
in order to estimate the systematic error related to the choice of the theoretical model for 
resonance production. In particular, the six models we use here are ALCOR [32], FRITIOF 
[28], Kaneta et al. [33], Letessier et al. [34], Stachel et al. [35] and UrQMD [36]. (See 
Appendix C for more details on the individual models.) We have also considered the AMPT 
2.11 A(mx)/A max simulation with string melting [37] that report on a non-thermal scenario 

9 



without an 77' mass modification. AMPT results show an interesting, although genuinely 
different dropping structure of A(mx) at low mx values. Comparison of this model to 
the data, as tersely overviewed in [9], indicate that AMPT cannot describe these data in 
a statistically acceptable manner, characterized by a y^/ndf = 102/13 corresponding to 
CL = 6.8 x 10~ 16 . We attribute the behavior seen in AMPT to a lower effective (ut) of the 
high mass halo resonances [9, 17]. 

In earlier simulations in [17], the resonance production was generated by an exponen- 
tial spectrum jV(mT-) = Ae~ mT ' Teff , the effective freeze-out temperature defined as T e ff = 
Tpo + m ( u t) 2 , with Tpo and (ut) being the freeze-out temperature and the average trans- 
verse flow, respectively. This has also been generalized, and a polynomial prefactor had been 
introduced in order to achieve a more realistic description of the direct production of reso- 
nances. We have fixed Tpo and (ur) to RHIC measurements [38]. Thus the distribution 
will follow the form of 

N(m T ) = Cm?e" mT/Te ", (10) 

where C is a normalization constant, and a — 1 — d/2, where d is the number of spatial 
dimensions of the expansion (hence 1 < d < 3 and a falls between —1/2 and 1/2) [19, 39]. 
The choice of a = corresponds to the case of Ref. [17]. 

In order to compare to RHIC data at mid-rapidity, we compute the effective intercept 
parameter of the 7r + -7r + correlations, A* (mi), with the definition of Eq. (4). All the contri- 
butions of Eq. (5) to the distribution were simulated in the mid-rapidity region, given 
by the pseudorapidity 5 cut of \y\ < 0.36 of the PHENIX acceptance. In the systematic 
checks we took into account that the STAR analysis [11] used a different pseudorapidity cut 
of \y\ < 0.5. 

For each model the fractions of the particles were computed by normalizing the multi- 
plicities to the total density 

Ptotal = Pcore + Phalo, Phalo = Pui + Prf + Pr) + Pk° ■ (H) 

The halo of pion production always had the same ingredients (from the decays of the u, 
r)', 1] and K° s ) 6 , while the core was composed of all the other resonances that were available 
In this work, pseudorapidity is denoted by y = 0.5 ln[( |p| +p z )/(|p| — Pz)]- In other manuscripts y is 
usually denoted by the variable 77 which in the current work is reserved for the pseudoscalar mesons rj and 
rf. 

6 Other long-lived resonances, such as the <j> meson, are checked to give a negligible contribution to the 
A* (tot)/ X" lax ratio, which translates to less thar|g% uncertainty to the m*, . 



in each particular model. The charged pion ttit spectra were obtained from a complete 
kinematic simulation of the decays above resonances using JETSET v7.4 [40]. The estimated 
systematic error arising from assigning each u decay product to the halo is given in Sec. V. 

The mechanism of the partial Ua(1) restoration implies that the 7/ would have a decreased 
effective mass in the hot and dense medium [17]. The number of the created rj' particles 
would then follow Eq. (10) with the modified mass and the freeze-out temperature respective 
to the rj' mesons, and the fraction of rj' mesons in the condensate is modeled with the 

N*, fml,\ a m v'- m > , 
N v > \m v > J 

formula, where the T con( i is the temperature of the condensate. 

As the escaping 7/ bosons are regaining their mass, they must lose momentum in order 
to fulfill the principle of energy conservation, i.e., 

* 2 , * 2 2 1 2 /io\ 

"V +Pt,t ) > = m v' +Pt, v ' ■ (13) 

(In the above equation the quantities with an asterisk denote the properties of the in-medium 
rf, while the ones without an asterisk refer to the free rf. According to the kinematical 
setup of both PHENIX and STAR measurements [10, 11], the longitudinal component of 
the rj' momentum is considered to be negligible here.) As a consequence, while rf bosons 
with ptp^ > \Jrn* r) , 2 — m v > 2 will follow the above distribution (with the effective mass m*, 
plugged in), the ones moving with a momentum less than this limit will be "trapped" in the 
medium until this medium is dissolved. Afterwards, rj'-s from the condensate are given a 
random transverse momentum, following Maxwell-Boltzmann statistics with a characteristic 
temperature B~ lm . 

Pl+P 2 y 

^K^^K 2 '"^ 1 ■ (14) 

with p 2 x + Py = Pt v >- Note that B is a systematically varied model parameter, and its best 
value is determined by the analysis of the \*{mi) data. 

V. THE 7/ MASS DROP RESULTS AND SYSTEM ATICS 

Out of the six parameters of the 77' spectrum and A*(t71t) simulation, the two most 
important ones are B~ l and m*,. They directly determine the shape of the observed dip 
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(depth and width, respectively). These parameters were considered as "fit variables," and a 
fine grid x 2 scan was carried out. The most probable values of the m*, and B~ l were looked 
for by the simulation of the A*(mx) spectrum and then determining the \ 2 from the fits 
to the data of the PHENIX and STAR measurements at each individual setup. The other 
four parameters were treated as "constants," although their variation played a role when 
determining systematic errors. 

As a default setup, the n + freeze-out temperature was taken from PHENIX inverse slope 
parameter fit 7 to be Tpo — 177.0 ± 1.2 MeV [38], and the conservative assumption of 
Tcond = Tpo was made. 8 The average transverse flow velocity was measured by PHENIX to 
be (ut) = 0.48 ± 0.07 (in relativistic units) [38]. This was cross-checked in our simulations 
with a 0.40 < (mt) < 0.60 scan with A (mt) = 0.01 steps at the best value of m*, and B -1 , 
and it was found that the most probable values were around (mt) = 0.50, in agreement 
with the PHENIX measurement. The polynomial exponent a = corresponds to a 2D 
expansion [19, 39], that had been proved to give a good description of data, used by PHENIX 
when obtaining (ut) and Tpo- However the a = +1/2 (ID) and a = —1/2 (3D) cases were 
also examined as systematic checks. 

The (-B -1 , mx) plane was mapped on a grid of 21 non-equidistant steps 9 for MeV < 
B- 1 < 350 MeV, and by 10 MeV equidistant steps 10 for the 10 MeV < m v > < 958 MeV 
region. A good simultaneous description of PHENIX and STAR data was obtained by 
certain values of B~ x and mx in the case of five out of the six models, as demonstrated in 
Figures 2-7. For each model, we show the confidence level (CL) map, the best result, the 
result without the mass modification, and the effect of m*, and B~ x scans around the best 
result. The standard deviation contours on (-B -1 , m^) plain for the best fits of A*(wt) / 'A™ aa: 
simulations to data are summarized for different models in Figs. 8, and 9 for a = and 
a = —1/2 respectively. 

As expected, a strongly reduced rj mass results in a dip in the low-mx part of the A*(?tit) 
distribution (and the smaller the mass, the deeper the dip is), while B~ x steers the width of 
this dip through smearing the momentum, while the (mt) parameter influences the overall 

7 Note that the values for Tpo and (ut) only include the statistical errors. The systematics were taken into 
account in our systematic studies that involved a wide range of different predictions on Tpo, 100 MeV < 
T FO < 177 MeV. 

8 The temperature of the condensate majorates the freeze-out temperature. As can be seen from our 

systematic checks, a freeze-out temperature that is lower than the condensate temperature will move 
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downward, i.e., toward "safety" for our conclusion. 

9 In fact, B^ 1 / 2 was mapped in 21 equidistant steps between and 600 MeV 1 / 2 . 
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FIG. 2. (Color online) Top left: A*(m T ) / '\f ax spectrum of the PHENIX and STAR data points, 
in comparison with Monte Carlo simulations based on the ALCOR model [32] at parameters from 
the best fit (B" 1 = 42 MeV, m*, = 490 MeV), for a = 0, T cond = 177 MeV, T FO = 177 MeV and 
(ut) = 0.48. Top right: Corresponding confidence level surface on the (-B -1 , my) plain. Only the 
region with acceptable fits (CL > 0.1%) is shown. Bottom left: A*(mT)/A^ ax spectra for different 
m*, values, -B -1 fixed to its best fit value. Bottom right: X^rwY) / \™ ax spectra for different B^ 1 
values, m*i fixed to its best fit value. In the latter two plots the PHENIX preliminary and the 
STAR Sinyukov measurements are also shown for comparison (as explained in Sec. III). 
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FIG. 3. (Color online) Results for the combined STAR+PHENIX data set with resonance multi- 
plicities from the FRITIOF model [28]. Explanation of the panels and other parameters are the 
same as in Fig. 2. Note that no values of B" 1 and m*, provided an acceptable fit, resulting in an 
empty confidence level plot. This model was therefore excluded from further studies. 

slope of the A t (ffiT) spectrum at higher tht- Such a role of (ut) was already observed in [17], 
and it is not detailed here. 

The systematic uncertainties resulting from the lack of precise knowledge of model con- 
stants were mapped by varying these constants within the range that is given by different 
measurements: Spectra with values of a = —1/2 and a = +1/2 were investigated as well 
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FIG. 4. (Color online) Results for the combined STAR+PHENIX data set with resonance multi- 
plicities from the model of Kaneta et al. [33]. Explanation of the panels and other parameters are 
the same as in Fig. 2. Best fit is at B~ l = 55 MeV, m*, = 530 MeV. 

as T F o = 100 MeV, T cond = 140 MeV and T cond = 220 MeV, representing the boundaries of 
the theoretically or experimentally acceptable region of these parameters. 

Besides the choice of the model parameters, there are systematic uncertainties resulting 
from uncertainties on A* (mr): 

• The STAR Edgeworth A* values were extracted from data of 0-5% centrality, compared 
to the 0-30% rapidity range of the PHENIX Sinyukov data set. We estimate that this 
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FIG. 5. (Color online) Results for the combined STAR+PHENIX data set with resonance multi- 
plicities from the model of Letessier et al. [34]. Explanation of the panels and other parameters 
are the same as in Fig. 2. Best fit is at B~ l = 86 MeV, m*, = 340 MeV. 

relative error does not exceed 9.8%. See details in Appendix B. 

• Although the u may give part of its contribution to the tail at the resolution available 
in STAR and PHENIX, we considered it as part of the halo. This introduces a relative 
systematic error of 7% on the measured A* [21]. 

• Moving the pseudorapidity cutoff from \y\ < 0.36 to \y\ < 0.50 adds another error, 
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FIG. 6. (Color online) Results for the combined STAR+PHENIX data set with resonance multi- 
plicities from the model of Stachel et al. [35]. Explanation of the panels and other parameters are 
the same as in Fig. 2. Best fit is at B" 1 = 86 MeV, m*, = 340 MeV. 

that is measured to be 3%. 

These errors are considered independent and are added up quadratically. The relative error 
on A*(mx)/A™ aa; corresponds to the same amount of relative error on m*, in the worst-case 
scenario. 11 

11 This relation was obtained from the STAR Edgeworth data set with the Kaneta et al. [33] multiplicities: 
when the lowest mx datapoint was moved by 10% first up and then down, the fitted m*, value changed 
by the same relative error. 
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FIG. 7. (Color online) Results for the combined STAR+PHENIX data set with resonance mul- 
tiplicities from the UrQMD model [36]. Explanation of the panels and other parameters are the 
same as in Fig. 2. Best fit is at B" 1 = 86 MeV, m*, = 400 MeV. 



The best fit values and their error determinations are summarized in Table I, while Ta- 
ble II shows the limits of acceptability, denned on the region, where CL > 0.1%. The 
combined PHENIX and STAR data cannot be described unless a significant rj' mass modi- 
fication is assumed, m*, < 680 MeV. Although this value includes the systematics from the 
model uncertainties, it is also subject to the uncertainties coming from the other sources 
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FIG. 8. (Color online) Standard deviation contours on the (B , mw) plain, obtained from 
\*(m T )/\™ ax of MC simulations for a = 0, T FO = T cond = 177 MeV and (u T ) = 0.48, based 
on different chemical freeze-out models, each fitted to the PHENIX and STAR combined data set. 
The framed band indicates the theoretically predicted range of 412 MeV-715 MeV [6], while the 
horizontal dashed-dotted line at 241 MeV indicates Weinberg's lower limit [42]. 
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FIG. 9. (Color online) Standard deviation contours on the (B , m^/) plain, a = —1/2. All the 
other settings are the same as in Fig. 8. 
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listed above. 12 On the basis of model simulations, studies and systematic checks, we find 
that a description of the combined PHENIX and STAR data set is possible with CL > 0.1% 
only if an in-medium 7/ mass modification of Am*, > 200 MeV is utilized. In other words, 
Am*, > 200 MeV at the 99.9% confidence level, corresponding to a more than 5-cr effect, in 
the considered broad model class. 

Note that the validity of our analysis relies on the correctness of published STAR and 
PHENIX data. These data, however, were not measured with the definite purpose to serve as 
a base for the search for the in-medium rf mass modification, where particular attention has 
to be payed to the momentum dependence of the particle identification purity and efficiency, 
especially at low-px regime. More detailed dedicated A*(mx) measurements together with 
additional analysis of the dilepton and photon decay channels of 7/ could help consolidate 
the findings reported here. 

VI. THE ENHANCED rf AND 77 SPECTRUM 

The dilepton spectrum has been measured recently in minimum bias Au+Au collisions 
a t \/ s nn = 200 GeV, and a large enhancement was observed in the low-invariant mass 
region m ee < 1 GeV [41]. Low transverse mass enhancement of the 77' and 77 production 
results in dilepton enhancement just in the considered kinematic range [6]. PHENIX recently 
reported a two-component transverse momentum spectrum in dilepton channel direct photon 
measurements [41]. The 77' and 77 spectra reconstructed here may serve as an input, e.g., for 
the simulations and evaluations of the dilepton spectra instead of the currently utilized 7/ 
spectra, based on m T scaling and hence not providing the possibility of taking into account 
an 77' mass reduction. 

In order to make a cross-check possible with measurements that are based on the non- 
pionic decay channels of the 77', we extract the transverse momentum spectra of the 77' and 
77 mesons from the Bose-Einstein correlation measurements at mid-rapidity in the reduced 
m*, scenario, for each of the successful resonance models. 

We compute the spectra two times for each model. First we use the original = 
958 MeV value. The obtained spectra clearly show the scaling. Then we use the 
12 We have determined these errors at the best mass fits, and applied the same absolute value in case of the 
mass limits. This is a conservative method since the higher the mass value is. to the lesser extent it is 
influenced by the same effect. 
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Acceptability boundaries of model fits 


Parameters 


Data set 


ALCOR 
[32] 


Kaneta 
et al. [33] 


Letessier 
et al. [34] 


Stachel 
et al. [35] 


UrQMD 
[36] 


a 


T CO nd 


Tfo 


m*, (MeV) 
S- 1 (MeV) 
CL (%) 


490+60 

42 
4.29 


530+5° 
55 
4.12 


340+50 

86 
6.35 


34Q+50 

86 
6.38 


400+5° 
86 
6.28 





177 


177 


X 2 /NDF 


1.83 


2.07 


1.72 


1.71 


1.72 








m* (MeV) 
S- 1 (MeV) 
CL (%) 


540+5° 

—40 

55 
2.80 


560+6° 

—30 

55 
3.35 


410+4° 

—40 

86 
6.07 


410+!° 

—40 

86 
5.97 


460+!° 

—40 

86 
6.14 


-0.5 


177 


177 


X 2 /NDF 


1.96 


2.07 


1.73 


1.73 


1.73 








m* (MeV) 




470 


210 












S- 1 (MeV) 
CL (%) 




55 
4.58 


86 
6.54 






0.5 


177 


177 


X 2 /NDF 




1.82 


1.71 












m* (MeV) 




620 


460 












S- 1 (MeV) 
CL (%) 




42 
2.26 


69 
5.86 









140 


177 


X 2 /NDF 




2.02 


1.74 












ml, (MeV) 




440 


200 












B- 1 (MeV) 
CL (%) 




69 
5.61 


104 
6.33 









220 


177 


X 2 /NDF 




1.75 


1.72 












m* (MeV) 




410 


240 












S- 1 (MeV) 
CL (%) 




145 
1.63 


145 
1.80 









177 


100 


x 2 M^ 




2.11 


2.09 













TABLE I. Fitted values of the modified rj' mass on the STAR+PHENIX combined data set, for 
different resonance models and parameters. The FRITIOF model has CL< 0.1% and therefore it is 
not shown here. The statistical errors are given by the 1-a boundaries of the fits, determined only 
for m*, and for the a = and a = —0.5 simulations. Best m*, and B -1 parameters for various 
systematic checks are shown in the last four rc^s. 



Data set 


ALCOR 
[32] 


Accepta 
FRITIOF 

[28] 


aility bour 
Kaneta 
et al. [33] 


Ldaries of i 
Letessier 
et al. [34] 


nodel fits 
Stachel 
et al. [35] 


UrQMD 
[36] 


Pe 

a 


irameti 

T C ond 


3rs 
Tfo 


PHENIX 

STAR 

combined 


0-750 
380-600 
380-590 


680-958 
none 
none 


0-720 
430-630 
420-620 


0-510 
190-450 
260-430 


0-500 
190-450 
260-430 


0-530 
260-500 
330-470 





177 


177 


PHENIX 

STAR 

combined 


30-770 
470-630 
450-620 


420-958 

none 
670-760 


50-730 
500-650 
490-640 


0-540 
300-500 
340-480 


0-540 
300-500 
340-480 


0-560 
360-540 
400-510 


-0.5 


177 


177 


PHENIX 

STAR 

combined 






0-690 
320-610 
340-590 


0-450 
0-390 
0-390 






0.5 


177 


177 


PHENIX 

STAR 

combined 






0-760 
560-690 
540-680 


0-450 
0-390 
0-360 









140 


177 


PHENIX 

STAR 
combined 






0-680 
270-580 
290-560 


0-410 
0-350 
100-320 









220 


177 


PHENIX 

STAR 
combined 






220-470 
360-470 
370-440 


30-310 

190-300 

200-280 









177 


100 



TABLE II. Acceptability boundaries of the modified rf mass on the PHENIX, STAR, and the 
combined PHENIX+STAR data sets, for different resonance models and parameters. A fit is 
considered acceptable if CL> 0.1%. FRITIOF fails completely on the STAR data set and also on 
the combined PHENIX+STAR data. All the other models require an m*, < 640 MeV excluding 
systematics. Acceptability boundaries for various systematic checks are shown in the last four 
rows. 

reduced rj' mass, and the corresponding B~ l value that provides the best description of the 
data in the frame of a given resonance model. These spectra break the mx scaling with an 
additional, steeper exponential-like part over the original rf spectrum, and produce most of 
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the rj' enhancement in the low-px region. The original "non-enhanced" spectra is normalized 
to the rj' multiplicity 13 predicted by Kaneta et al. [33], and then the enhanced spectra is 
scaled relative to it. Left panel of Fig. 10 shows both the original, m-p-scaling spectra, and 
the enhanced version with m*, < for the resonance ratios of Ref . [33] . A comparison of 
enhanced rj' spectra of all resonance models is given on the right panel of Fig. 10, and the 
fitted spectrum parameters are listed in Table III together with the rj' enhancement factors 
for each particular model of resonance ratios. 




FIG. 10. (Color online) Left: Reconstructed tot spectrum of the rj' mesons using the resonance 
ratios of Ref. [33]. Lower (blue) part indicates the scenario without in-medium rj mass reduction, 
upper (red) part the enhancement required to describe the dip in the low tot region of A*. Right: 
Comparison of reconstructed tot spectra of the rj' mesons from different models. All spectra are 
normalized to the rj' multiplicity as given by the model of Kaneta and Xu [33] , and then fitted with 
a double exponential. The shaded (yellow) band represents the total error. The fit parameters 
are summarized in Table III. Above tot — m,qi = 1 GeV, all models result in very similar values, 
corresponding to an approximate mx-scaling. This figure indicates that the violation of this tot- 
scaling is model dependent, and suggests that dilepton measurements may provide additional 
constraints for the model builders. 

13 The normalization of the original spectra is performed the following way: First the numerical integral 
of the unmodified -^-^^ distribution is computed. Then the average number of rj' mesons per event is 
divided by this integral, and the histogram is scaled by the resulting number. 
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Resonance model 


a 


b 


c 


d 


fr}' 


No enhancement 


0.82 


2.72 


(fixed) 


(fixed) 




ALCOR [32] 


2.30 


2.98 


62.4 


23.5 


43.4 


Kaneta et al. [331 

L J 


2.21 


2.94 


32.4 


18.7 


25.6 


Letessier et al. [34] 


2.91 


3.14 


80.1 


12.8 


67.6 


Stachel et al. [35] 


2.85 


3.13 


80.0 


12.8 


67.6 


UrQMD [36] 


2.75 


3.07 


52.5 


12.7 


45.0 



TABLE III. T]' enhancement (Jy) and fit parameters of the spectra for different models of resonance 
abundancies. The spectra are obtained using the most probable B" 1 , m*, parameters, and then 
fitted with a double exponential function y = ae _b '- mT_ "V) + ce ~ d ( mT ~ m v'\ The spectrum without 
enhancement is computed with the Kaneta et al. [33] resonance ratios, with m*, =m^i . 

Since the spectrum of the rj mesons can be directly compared to measured data, it also 

serves as a consistency check. A comparison of enhanced rj spectra of all resonance models 

is given on Fig. 11 (left), and the fitted spectrum parameters are listed in Table IV. A 

connection between the rj' enhancement and the 77 enhancement can be expressed as 

N* ( N* \ N 1 

^tr l+ \t- l )t m ^ v+ ^ ] (15) 

with the last term, the total rj to rj branching ratio being approximately 65.7%. For models 
[32-35] we have also plotted the absolutely normalized rj spectra in Fig. 11 (right) com- 
pared to the measured PHENIX rj spectra. It is obvious from the plot that the computed 
enhancement affects only the < 1 GeV part of the rj spectrum. 

Considering all errors, 6.01 < N r) > < 258 and 1.56 < N v < 15.4. 14 Estimations using 
these enhancement factors indicate that the observed in-medium rj' mass drop is indeed a 
promising candidate [9] to explain the dilepton excess reported by PHENIX in Ref. [15]. 
14 The upper bounds on N v > and N n were calculated the following way: The (-B -1 , m*,) values of all above 
mentioned setups of input model, a, Tpo, T con d were considered along the 1-a contour, and then the 
m*i value was shifted upward with the corresponding non-model-dependent systematic error. Then the 
•q' (or rj) spectra were plotted for these (B , m*,) pairs, and the one corresponding to the maximum 
enhancement was considered. The lower bound was computed similarly, with the m*, values shifted 
downward and the spectra with the minimal enhancement taken. The maximal enhancement is given by 
Ref. [34], while the minimal is by Ref. [33]. (Ref. [28] is not considered as FRITIOF fails to describe the 
STAR and the combined PHENIX and STAR data sets.) 




FIG. 11. (Color online) Comparison of reconstructed uit spectra of the r\ mesons from different 
models. Left: All spectra are normalized to the 7/ multiplicity as given by the model of Kaneta and 
Xu [33], and then fitted with a double exponential. The shaded (yellow) band represents the total 
error. The fit parameters are summarized in Table IV. Above mx — = 1 GeV, all models result 
in very similar values, corresponding to an approximate mi-scaling. Right: Absolute normalized 
spectra from input models [32-35] are compared to PHENIX 200 GeV central Au+Au collision 
measurements [43]. 



Resonance model 


a 


b 


c 


d 


fv 


No enhancement 


14.6 


3.38 


(fixed) 


(fixed) 




ALCOR [32] 


14.6 


3.40 


97.0 


17.8 


5.25 


Kaneta et al. [33] 


14.6 


3.38 


54.9 


16.2 


3.47 


Letessier et al. [34] 


14.6 


3.38 


84.1 


16.9 


4.75 


Stachel et al. [35] 


14.5 


3.38 


89.2 


17.0 


4.97 


UrQMD [36] 


14.6 


3.41 


148 


17.9 


7.49 



TABLE IV. rj enhancement (f^) and fit parameters of the spectra for different models of resonance 
abundancies. The spectra are obtained using the most probable -B" 1 , m*, parameters, and then 
fitted with a double exponential function y = ae ~ b(mT ~ r 'V) + C e~^ mT " m v). The spectrum without 
enhancement is computed with the Kaneta et al. [33] resonance ratios, with m*, =m^i . 



25 



VII. CONCLUSION 



The best simultaneous description of STAR and PHENIX HBT data is achieved with 
an rf mass that is dramatically reduced from 958 MeV to ± 42 MeV in the 

medium created in central Au+Au collisions at RHIC. The first error here is the statistical 
one from the fit, the second error is from the model and parameter choices, while the third 
is the systematics from the centrality selection, the resolvability of u decay products and 
the pseudo-rapidity cutoff. Note that the dominant error corresponds to the selection of the 
model for the hadronic multiplicities. 

Note that the best estimated value for the modified 7/ mass does not differ significantly 
from the range of 412 MeV — 715 MeV, predicted by the quark model calculations for the 
Ua(1) symmetry restoration [6]. In fact, it is slightly below this range, but above the lower 
limit of y/3m n by Weinberg [42]. Hence the mass reduction effect seems to be already at 
maximum at ^s NN = 200 GeV central Au+Au collisions. As a consequence, one may expect 
a saturation of this effect if the bombarding energy is further increased up to LHC energies 
of y/s NN = 10 TeV. In Appendix B we also noted an interesting centrality dependence of the 
order of 9.8%, that suggests that the in-medium rf mass decrease is slightly larger in more 
central collisions and suggests that more detailed centrality and transverse mass dependent 
measurements of the Edgeworth and other model-independent estimates of the extrapolated 
intercept parameter A* are necessary to pin down this effect, which at present is part of the 
above given 42 MeV systematic error. 

Not only did we investigate the best value for the in-medium mass of 7/, but also the 
minimum mass modification that is required to describe the data. Based on the combined 
STAR and PHENIX data, and from the systematic investigation of various resonance mul- 
tiplicities and model parameters, we conclude that in y/s NN = 200 GeV central Au+Au 
reactions the mass of the 77' meson is reduced by more than 200 MeV, at the 99.9 % confi- 
dence level in the considered model class. This result was briefly summarized and published 
in Ref. [9]. A similar analysis of NA44 S+Pb data at top CERN SPS energies provided 
no evidence of an in-medium rf mass modification [17]. For further details on the extrap- 
olation techniques and on a summary of earlier results on a correlation search for partial 
Ua(1) symmetry restoration we recommend the review paper [44]. Theoretical results and 
earlier experimental searches for in-medium mass modifications of hadrons were summarized 
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recently in Ref. [45]. 

Our positive results on a significant, indirectly observed in-medium rj mass modification 
should revitalize theoretical interest in other signatures of partial Ua(1) and chiral symmetry 
restoration in heavy ion reactions and also should be cross-checked against other observables 
like the enhancement of low-mass dileptons and photons in the same reactions. More detailed 
and more precise experimental data are needed on the intercept parameter of Bose-Einstein 
correlations of pions at low at various bombarding energies, system sizes and centralities 
to fully understand the nature of partial Ua(1) symmetry restoration. Detailed analyses of 
other decay channels of the if and t] mesons, e.g. dilepton measurements are required to 
confirm our observations on the onset of this effect. 
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Appendix A: STAR Edgeworth A* 

STAR used a sixth order Edgeworth expansion with the even order terms only, assuming 
that the source is a symmetric and analytic function. The STAR measurements for the A^ 
and k values are given in Table V. Plotted values of the Gauss A and the fourth and sixth 
order Xe, compared to our A* calculations, are shown in Fig. 12. In this Appendix we outline 
the computation of A*, using the notation of STAR [11]. Although the Hermite polynomials 
were defined differently from Eq. (8) of Ref. [11], the difference cancels from the equations 
quoted below. The o, s, I indices denote the out, side and long directions as defined at the 
same place. Then the A* is computed as 

A* =\e (l + \ K o,A — igKofi) (l + |k Si 4 — ^^Ifi) (l + §«7,4 — i§ K s,6) , (Al) 
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and its error is given as 
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9ki,6 ' J \dKi >6 ' J \d\ E 
The result of the numerical computations are summarized in Table V. 
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FIG. 12. (Color online) HBT Gaussian A, fourth and sixth order Edgeworth A^; parameters from 
Ref. [11], and the computed A* parameters are shown for the four transverse momentum bins. 
(Note that, for the sake of visibility, the points are slightly shifted horizontally to the left or right.) 



Appendix B: System size, energy and centrality dependence 

The STAR Edgeworth A* (m-r) data are given for the 0-5% most central data, while 
PHENIX carried out the pion correlation analysis on the 0-30% centrality range. The 
behavior of the A* (%) curve is not necessarily the same for different centrality classes, 
resulting in a systematic error in our analysis. Fortunately, STAR provided a centrality de- 
pendent Gaussian A measurement, that we could use to estimate this error. The A(mx) / \ ma x 
values for the centrality classes between 0% and 30% are summarized in Fig. 13. The relative 
error is determined to be less than 9.8 % given by the the difference of the first A(m T )/A max 
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Pt 
(MeV/c) 


150-250 


250-350 


350-450 


450-600 


A 


0.30+0.01 


0.42±0.01 


0.45±0.01 


0.47+0.01 


\e 


0.23+0.01 


0.35±0.01 


0.41±0.01 


0.44+0.01 


Kfy A 
J - * 


0.53±0.11 


0.45±0.10 


0.20±0.11 


0.22+0.13 




0.83±0.39 


0.53±0.38 


0.63±0.44 


-0.84+0.53 


Kg 4 


0.99±0.10 


0.79±0.10 


0.16±0.11 


-0.07±0.13 


Ks,6 


3.07±0.35 


3.21±0.37 


1.71±0.44 


1.80±0.51 


Kl,4 


1.32±0.07 


0.70±0.07 


0.54±0.09 


0.32±0.11 




-1.76±0.29 


-2.82±0.29 


-2.41+0.35 


-2.12±0.43 


A* 


0.307±0.015 


0.433±0.015 


0.456+0.016 


0.475+0.018 



TABLE V. HBT Gaussian A parameters, Xe and K%^ n fit parameters of the sixth order even Edge- 
worth expansion for the 5% most central events as found by STAR [11], used to obtain the also 
shown A* parameters. (We kept one extra decimal here in order to avoid fitting errors later due to 
the coarse roundoff of the uncertainties.) 

datapoints of the 0-5% with respect to the 20-30% centrality classes. This is a conservative 
estimation of the difference between 0-5% and 0-30% centrality data. 

According to this centrality and transverse mass dependent Gaussian STAR data on 
the intercept parameter, the depth of the low-mx dip deepens in the case of more central 
collisions. This suggests that the restoration of the Ua(1) symmetry is more complete in 
case of more central collisions. Therefore it may be useful to measure also the more relevant 
Edgeworth A* (m-p) data in different centrality classes too. 

System size and energy dependence is shown in Fig. 14, where the A(/cT)/A maa . values 
are compared for NA44 S+Pb collisions, and STAR Au+Au and Cu+Cu collisions. The 
plot indicates that the mass modification effect seems to be maximal in 200 GeV Au+Au 
collisions, followed by 62 GeV Au+Au collisions, 200 GeV Cu+Cu and 62 GeV Cu+Cu 
in that order. Although in each of the STAR cases a positive signal is apparent, we also 
observe the lack of saturation of A(/cT)/A maa; at higher momenta, which we attribute to the 
changing (decreasing) amount of the radial flow with decreasing energy and system size. 
Such a dependence on the amount of the radial flow has been first pointed out in Ref. [17]. 

29 



XIX max from J. Adams et al., STAR Coll., Phys. Rev. C71, 044906 (2005) 
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FIG. 13. (Color online) STAR HBT Gaussian A(m T )/A max values for different centrality classes 
between 0% and 30% [11]. The points with different centralities were measured in the same wit 
intervals, but are slightly shifted to left and right on this Figure for a better visibility. Note 
the trend that more central data correspond to a larger hole in the lowest transverse mass bin, 
suggesting a slightly larger in-medium rj 1 mass decrease in more central collisions. This effect is 
part of the systematic errors given in the conclusions. 

Appendix C: Models 

In this section we provide a brief description — both general and analysis specific — of the 
six different models used to simulate the number of each important resonance that decays 
into pions, as well as the number of primordial pions. The fractions of pions from different 
sources are compared to the total number of particles in Table VI. These values are fed 
into our simulations as an input. Note that not all core resonances are listed. With the 
assumption of a fixed the A»(mx) is solely determined by the number of the halo 

resonances (u, rj, rf, K° s ), and by the total number of pions. The primordial pions can 
replace those pions that come from fast decays without any effect on A»(mx). Generally, if 
a model contains a large number of exotic short lived resonances, it predicts less primordial 
pions. 
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FIG. 14. (Color online) STAR Gaussian \{k^)/\ max values for RHIC collisions of different syst 
and energies [48], compared to SPS NA44 measurements [49]. 







Fractions from model 




Particles 


ALCOR 


FRITIOF 


Kaneta 


Letessier 


Stachel 


UrQMD 




[32] 


[28] 


et al. [33] 


et al. [34] 


et al. [35] 


[36] 


primordial tt + 


0.4910 


0.2095 


0.7396 


0.3059 


0.3333 


0.2395 


P 


0.1100 


0.3058 


0.0651 


0.0259 


0.0370 


0.0045 


A 




0.0846 


0.0088 


0.0080 




0.0069 


short lived K* 




0.1351 


0.0124 


0.0143 




0.0014 


£ 




0.0153 


0.0040 


0.0017 




0.0066 


S* 




0.0098 


0.0054 






0.0016 




0.1104 


0.1023 


0.0209 


0.0233 


0.0296 


0.0073 


V 

long lived 

rf 


0.0453 
0.0067 


0.0516 
0.0577 


0.0602 
0.0089 


0.0383 
0.0032 


0.0360 
0.0031 


0.0233 
0.0050 




0.0717 


0.0283 


0.0746 


0.0601 


0.0513 


0.0287 



TABLE VI. Resonance ratios from different models. 
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ALCOR [32]: Hadron multiplicities — especially for strange particles — are calculated in 
the framework of this algebraic coalescence rehadronization model, which counts for redis- 
tribution of quarks into hadrons for relativistic heavy-ion collisions. Resonance ratios were 
taken from recent calculations for the RHIC collisions [46]. Since the mass eigenstates of the 
strange-antistrange mesons are not directly pinned down, the Kaneta-Xu [33] predictions 
were used to fix the N v r /N v ratio for the ALCOR model. 

FRITIOF [28] is a Monte Carlo program that implements the Lund string model for 
hadron- hadron, hadron-nucleus and nucleus-nucleus collisions. Resonance ratios were com- 
puted with the FRITIOF model by simulating of 1000 events using RHIC central Au+Au 
y/s NN = 200 GeV parameters. FRITIOF was excluded from examination when drawing the 
conclusions of this analysis, since it was unable to describe the STAR A*(mx) / \™ ax data set, 
nor the combined STAR + PHENIX data set, to a CL > 0.1% at any setup of the model 
parameters. 

Kaneta-Xu [33]: In the case of this model, Eq. (1). of [33] was used the following way: 



where m 8 is the mass, Jj is the spin of the particle. K 2 (x) denotes the second order modified 
Bessel function. For the sake of simplicity the temperature of the chemical equilibrium 
T c h = 160 MeV, the strangeness saturation factor 7 s = 1, the quark potentials \x q = 10 MeV 
and \i q — were considered, all in consistency with the published PHENIX and STAR 
measurements in y/s NN = 200 GeV Au+Au collisions [47]. (The relative minuteness of \i q 
guarantees that the charge factor is negligible for even those rare particles, where it is not 
exactly 0). 

Letessier-Rafelski [34]: The model of these two authors studied soft hadron production 
in relativistic heavy ion collisions in a wide range of reaction energy, 4.8 GeV < y/s NN < 
200 GeV, and made predictions about yields of particles using the statistical hadronization 
model. Particle yields of the central events were taken from Table 4 of [34]. 

Stachel et al. [35]: This statistical "fireball" model treats the system as a grand canon- 
ical ensemble with the temperature and the baryon chemical potential as free parameters. 
Particle yields of the central events were taken from Sec. 2.2. Table 1 of [35]. 

UrQMD [36]: The ultra-relativistic quantum molecular dynamics model is a microscopic 
model used to simulate (ultra)relativistic heavy ion collisons in the energy range from Bevalac 




(CI) 
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and Heavy Ion Synchrotron (SIS) up to Alternating Gradient Synchrotron (AGS), SPS and 
RHIC. Resonance ratios were computed with the UrQMD model by simulating 1000 events 
using RHIC central Au+Au y/s NN = 200 GeV parameters. 
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